General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



N76- 309?0 


' lAb/ ION '■o-r’-'FAT'J?»;S IS HIP-1 
I) Fpr>v .:hai £.-5<ch;sgi 

”‘ICAL !'r;S5lOS «r> XCI’^A (SuSA) 41 p HC 
14.00 JSCL 201 Unci IS 

03/73 5U4lb 


NASA TECHN ICAL 

MEMORANDUM 


NASA TM X-73471 


rs» 

I 

X 



< 

< 


ION ItMPERATURES IN HIP-1 AND SUMMA FROM 
CHARGE -EXCHANGE NEUTRAL OPTICAL EMISSION SPECTRA 


by R. W. Patch and M. R. Lauver 
Lewis Research Center 
Cleveland, C -‘.4' 35 


TECHNICAL PAPER presented at 

Conference on Diagnostics of High Temperature Plasmas 
sponsored by the Division of Plasma Physics 
of the American Physical Society 
Knoxville, Tennessee, January 7-9, 1976 



i:-a6oi» 


ION TIiMPERATURiiS JN HIP-1 AND SUMMA FROM 
CHARGE-EXCHANGE NEUTRAL OPTICAL 
EMISSION SPECTRA 
by R. Mf. Patch and M. R. Lauver 
Levis Research Center 


SUMMARY 

Ion temperatures were obtained from observations of the H^ , D , and 
He 587.6 nm lines emitted from hydrogen, deuterium, and helium plasmas in 
the SUMMA and HIP-1 mirror devices at Lewis Research Center. Steady-state 

c * 5 discharges were formed by applying a radially inward dc electric 
field between cylindrical or annular anodes and hollow cathodes located 
at the peaks of the mirrors. The ion temperatures were found from the 
Doppler broadening of the charge-exchange components of spectral lines. 

A statistical method was developed for obtaining scaling relations 
of ion temperature as a function of current, voltage, and magnetic flux 
density . The method was applied to hydrogen with ion temperatures between 
169 and 85 ^ eV, deuterium with temperatures between 266 and 65 O eV, and 
helium with temperatures between 597 and 1*927 eV. 

Derivations are given that take into account for the first time tri- 
angular monochromator slit fur.ctions, loss cones, and superimposed charge- 
exchange processes. In addition, the Doppler broadening was found to be 
sensitive to the influence o. drift on charge-exchange cross section. 

The effects of finite ion-cyclotron radius, cascading, and delayed emission 
ore reviewed. Including the effect of slit function and loss cones re- 
duced the ion temperatuire. In the zero-ion-cyclotron-radius approximation 
including drift increased the ion temperatui*e. In the finite-ion-cyclotron- 
radlus approximation an additional physical phenomenon tended to decrease the 
ion temperature at high ion temperatures. 


INTRODUCTION 

This report presents results of research undertaken to measure and 
obtain scaling relations for ion temperatures in the Hot Ion Plasma 1 
(HIP-l) apparatus and the Superconducting Magnetic Mirror Apparatus 
( SUMMA) at NASA Lewis Research Center. Both HIP-1 and SUMMA are magnetic 
mirror facilities. HIP-1 has copper magnet windings whereas SUMMA has 
superconducting magnet windings. In both facilities a steady state dis- 
charge is formed by applying a radially inward dc ilectric field between 
cylindrical anodes and hollow cathodes located at the peaks of the mirror. 
The perpendicular electric and magnetic fields cause the electrons and 
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Ions of the plasma to drift azlmuthally with approximately the same 
velocity, corresponding to an ion drift energy of as mtich as several 
kilovolts. In addition, it appeal’s that turbulence occurs, which results 
in Ion heating (ref. l). Die ion temperatures in such devices are many 
times the electron temperatures. 

The approach used in measuring ion temperatui'es was to observe emis- 
sion from the Doppler-broadened charge-exchange neutral spectra of atomic 
lines. The method used in this report for obtaining ion temperatures from 
the measured width (or second moment) of such lines was given in refer- 
ence 2 but was extended to deuterium for the present work. 

Statistical methods were used to obtain improved scaling relations 
for ion temperature from and He 587.6 nm SUMMA data (ref. 2) and for 

similar data from H and D HIl'-l observations. The SUMMA data were ob- 

a a 

tained for uncooled tungsten hollow cylindrical cathodes and cooled anodes, 
and the HIP-1 data were obtained for tungsten-coated water-cooled copper 
cathodes and ring anodes. 

The scope of the scaling relation study is limited to average ion 
temperature in pure D 2 , and He plasmas at the magnetic midplane. The 

averages are obtained by observing the emission spectra across a diameter 
and are Inherently heavily weightea in favor of the center of the plasma, 
where most of the light intensity originates. No actempts were made to 
invert the observed spectral intensities to get intensity and hence ion 
temperature as a function of radius because an Abel inversion is not ap- 
plicable due to the anisotropy of the spectral emission. 

Also contained in this report are some advances in the theory of 
broadening of charge-exchange neutral spectra from rotating plasmas in 
mirror machines. In the future these advances will be used to obtain 
more accurate ion temperatures than can be obtained by the method in 
reference 2. 

Closely related experimental efforts on charge-exchange neutral 
Doppler broadening were presented in references 3 to 5. An influence of 
finite lifetime of excited states was observed in the experiment of ref- 
erence 6. A plasma model which includes configuration as well as velocity 
space (ref. 7) has been recently adapted to Doppler broadening. The in- 
fluence of lifetime of excited states, finite cyclotron radius, and azi- 
muthal drift was studied in reference 8. 

The influence of drift on the ov correction to line shape is pre- 
sented herein as an extension to reference 2. The influence on broadening 
of cut-off of velocity distribution due to loss cones is also discussed. 

The authors are indebted to Dr. Norman To Ik of Bell Telephone Labor- 
atories for measuring the r itio of excitation cross sections of the 
He 587.6 and He 667.8 nanometer lines. 
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EXPERIMEHTAL APPARATUS 
SUI>!MA 

Ma«^uetB . - The Superconducting Magnetic Mirror Apparatus (SUI<C4A) 

(ref. 10) includes four superconducting solenoldal magnets of which the 
inboard two were used for present study (see fig. l). Each magnet has 
an inner, a middle, and an outer winding, which were Independently pow- 
ered. The two inboard m«ignets alone produce maximum steady flax densities 
of about 5-2 T and 3.3 T at the mirrors and midplane, respectively. The 
average mirror ratio R was 1.33 (symbols are given in appendix A). 

Plasma test section . - A schematic view of the plasma test section 
is shown in figure 1. The discharge chaxi.ber is a cylinder 3.73 m in 
length and 36.6 cm in diameter made from type 30U stainless steel. Two 
23. ^ cm diffusion pumps provided a base pressure of 5*10"^ torr. 

Electrode assembly . - The electrode arrangement is also shown sche- 
matically ir; figure 1. More detail is shown in figure 2. Gas is intro- 
duced through the uncooled tungsten cathodes (1.27 cm o.d. ) to produce 
a hollow- cathode discharge. The water-cooled copper anode rings have 
2.9 cm and 3.1 cm inside and outside diameters, resepectively. They are 
electrically grounded to the test section walls. The electrically 
floating shield protects a cathode holder from heavy particle bombard- 
ment and helps inhibit arcing from the cathode to the electrically 
grounded back flange. The shield was constructed of tungsten-coated 
copper and cooled by water. The cathodes are connected in parallel to 
a single dc power supply. 


HI P-1 

Magnets . - In this apparatus two water-cooled copper solenoidal 
magnets are used to produce a maximum steady flux density of l.l6 T at 
the midplane and 2.13 T at the mirrors for a fixed mirror ratio R of 
1 . 82 . 

Plasma test section . - The discharge chamber is a cylinder 1.63 m 
in length and 25 cm in diameter made from type 30U stainless steel. 

Two 33.6-cm and two 15.2-cm diffusion pumps provide a base pressure of 

-7 

5x10 torr. 

Electrode assembly . - The two electrode assemblies, of which one 
is shown in figure 3, are placed near the mirrors in an arrsuigement 
similar to SUMMA. The gas supply to the hollow cathodes is similar to 
SUMMA, but the ma+erials, cooling and dimensions of the electrodes differ. 
The cathodes (W.C. no. 3 in ref. 9) have inside and outside diameters 
of 2.3^ cm and 14.1*5 cj., respectively. The auioies have inside and outside 
diameters of 6.35 cm and 0.25 cm, respectively. The cathodes, anodes, 
and floating shields are all water-cooled copper with plasma flame- 
sprayed coatings of tungsten (see fig. 3). Electrical connections are 
similar to 3UMMA. 
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Di agnostic Equipment 

The same apparatus for emission spectroscopy was used In SUMMA and 
HIP-1 and Is shown In figure k. It was located at the test-section mid- 
plane and viewed the plasma perpendicular to the magnet axis. The ver- 
tical centerpleuie (a vertical plane contoinin- the magnetic socls) was 
focussed on the entrance slit of the monochromator with a lens. Between 
the lens a.d the entrance slit the beam was rotated 90® with a beam rotator. 
A stop was provided to reduce the vertical height of the beam. The beam 
cross-section at the vertical centerplane was detemined by the size and 
shape of the entrance slit. The beam cross-section was approximately 
2 mm in the vertical direction and 2.6 cm parallel to the magnetic axis. 

An adjustable screw at the base of the monochromator table was used to 
make vertical scans. For the work reported here the monochromator was 
focused across the centerline of the plasma so the quantity y in fig- 
ure 1* was zero. The grating monochromator employed an f/8.6 1/2 meter 
Ebert mounting with curved slits. The reciprocal linear dispersion was 
1.6 nm/mra in first order. A photonultipller with extended S-20 photocathode 
was magnetically shielded and cooled with air that came from a vortex tube 
refrigerator. 

In SUMMA a floating electrical probe equipped for rapid traverse of 
the plasma was Installed for several runs. 


TEMPERATURE SCALING REUTIONS 
Theory 

Scaling relations of the form 


T^ . 10 ° I ^3 3 

were found in HIP-1 (ref. 9), where T^ is ion temperature at the magnetic 

midplane, I is total current into the cathodes, V is electrode voltage, 

B is mldplane magnetic flux density, and b^ through b^ are constants. 

The form of equation (l) was found by the process of holding two independent 
variables fixed and varying the third. If equation (l) is inexact, such a 
procedure may bias the data. Use of a regression analysis makes the pro- 
cedure of reference 9 unnecessary and makes it possible to take the data 
without biasing it and without the experimental complication of holding 
two independent variables fixed. 

To do a regression analysis, we write equation (l) in logarithm form: 
log T^ = Bq + b^ log I ^2 V + b^ log B (2) 

This is the standard form of the regression equation for triple linear 
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correlation (ref. U), where log I, log V, and log B are Independent 
variables, through b^ are coefficients of net regression, and 

log Tj Is the dependent variable. The coefficients of net regression 

may be found by tne methods in reference 11. 

The use of the above method makes it possible to obtain scaling 
relations from old data where the independent variables were not held 
fixed. 

To estlrAte the reliability of the regression coefficients and the 
resulting values of T^ , a statistical analysis (refs. 11, 12, 13) is 

useful. The multiple regression euid statistical analyses were performed 
with three programs for a high-speed digital computer. Only one of these 
has been documented (ref. 12). 


Results and Discussion 


Hydrogen and helium were run in SUI«IMA. In HIP-1 these two gases 
and deuterium were run. Test conditions are given in table I. The 
data for HIP-1 are for only one anode ring installed. Observations were 
made of the H^, D^, and He 587.6 nm lines. A typical observation of 

H^ is shown in figure 5 from SUMMA. No forbidden satellites due to 

plasma turbulence were identified. 

The runs differed in the following particulars from the HIP-1 runs 
reported in reference 9: l) In this report gas flow for both apparatuses 

was through both cathodes whereas in reference 9 it was through one 
cathode only, 2) In this report in order for the statistics to be r »i ' g- 
ful each data point represents one set of experimental readings wher.«-; 
in reference 9 each data point represented an average of several sets of 
experimental readings, 3) In this report no attempt was made to take data 
at any particular point if the anode current varied cyclically, whereas 
in reference 9 all data used in detemining scaling relations was taken 
at the lowest auiode current in the cycle, k) In this report no attempt 
was made to hold two variables constant and vary a third, but this pro- 
cedure was followed in reference 9. 

In converting Doppler broadening tc ion temperature, the methods 
of reference 2 were used because numerical methods for handling large 
amounts of data have not been developed as yet for the case involving 
drift and loss cones. It was assumed that the only significant reactions 
were 

H* + H 2 ^ H* + H^ (I) 


d"*^ + D« 
f 2 


* + 
Df .D2 


(II) 
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HeJ ♦He - He* ♦ He"^ (ill) 

where subscript f stands for fast and superscript * stands for elec> 
tronlcally excited. The only corrections made to the second moments of 
the charge-exchange components of the spectral lines were for fine stjruc- 
ture, magnetic splitting, and slit function (ref. 2) and the total cor- 
rections never exceeded 11 percent. Ion temperatures were then found 
from figure 6, which is the same as table I of reference 2 except/ that 
D and D. have been added. 

On a 

♦ ♦ 

The scaling relations for H and H^ in SUMMA are given in fig- 
ures 7 and 6. The scaling relation for H in HlP-1 is given in figure 9. 

For D'*’ in HIP-1 the voltage was held constant, so only a double correla- 
tion W 61 S possible. It is given in figure 10. For He*^ in HIP-1 there 
was too small a sample to obtain a reliable scaling relation c'' any kind. 
The scatter in figures 9 and 10 is greater than in reference 9 because of 
differences 2) and 3) mentioned above. 

A statistical summary of the T^ scaling relations is given in 

Table II. For H* in &UMMA and HIP-1 the correlations with I, V, and B 
were all highly significant (definitions of "significant" and "highly 
significant" are given in ref. 13). For He'*' in SUMMA the correlation 
with V was highly significant, but the correlations with I and B 
were not significant. This probably means that a larger sample is re- 
quired, or I and 6 must be varied over wider ranges, or experimental 
variables must be measured more accurately. For HIP-1 with D* the 
correlation with I was highly significant whereas the correlation with 
B was only significant. 

Standard deviations • of coefficients of net regression are also given 
in Ted>le II. With these the scaling relations may be written: 


„ ^0.52t0.10,,1.61*±0.fJ5,,-1.30±0.22 

% 1 V D 

(SUMMA H"^) 

(3) 

„ T-0.11t0.2U,,1.36t0.05o-0.29t0.22 

^ 1 V Jj 

(SUMMA He"^) 

(M 

^ jO. 33i0. OUyl. 50t0. lOg-0. U9t0. 10 

(HIP-1 h'*’) 

(5) 

^ j0.32t0.053^0.30t0.12 

1 D^) 

(6) 


Equation (5) does not agree exactly with reference 9, where for 
slightly different conditions it was found 

Tj^ -V (hip-1 h'*') 


(7) 
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The difference in I and V exponent! is within the experimental error. 

The cause of the discrepancy in B exponents has not been investigated. 

The last column of Table II contains the standard error of estimate 

S of log.. T. . The relative standaurd error of estimats can be found 
z ”10 1 

from the last column and the relation 


which gives 



S In 10 « 2.30 S 
z z 


0.0685 (SUMMA H'*’) 


O.lOU (SUI-WA He*) 


0.112 (HIP-1 H*) 


0.136 (HIP-1 D*) 


( 6 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 


Hence for H* in SUMMA the ion temperatures deviated on the average about 
6 . 85 ^ from the value estimated from the scaling relation (together with 

» 0 >- 

In order to compare H data from HIP-1 and SUMMA, the SUMMA data 
were replotted in figure 9 even though the correlating parameter is not 
optimum for the SUMMA data. It can be seen that the uncooled tungsten 
cathodes in SUMMA gave higher temperatures than the cooled tungsten- 
coated copper cathodes in HIP-1. This is partly due to the difference 
in cathode diameters as was verified by running cathodes of smaller diam- 
eter in HIP-1 (ref. 9). It may also be partly due to a higher neutral 
density in HIP-1 causing higher charge-exchange power loss (see ref. 9). 
However, UIP-1 plasmas had leurger diameters than in SUMMA, which is 
desirable. 

The results in figures 7, 9, and 10 should be viewed with some 
caution, as order-of-magnltude nonequilibrium composition calculations 
indicated that there could be 15 percent as much B* as H* , so the 


reactions 
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^2f * ^2 * ^ 

and •imllar reactions for may be significant although they were 

neglected. However, because the shape of the cross section curve for 
reaction IV resembles the shape for reaction I, the resulting error in 
should not be serious. 

The line shape of from SUMMA agrees well with the line shape 

calculated from equation (B7) of reference 2 as shown in figure 11. For 
these conditions the fine-structure , magnetic-splitting, and slit- 
function corrections were negligible. The effect of azimuthal drift was 
neglected in computing the curve labeled "calculated" in figure 11. 


IMPROVED THEORY OF LINE SHAPES AND SECOND MOMENTS 

A theory of line shapes and second moments of charge-exchange 
neutral spectra was given in reference 2. In this section the treatment 
of reference 2 ij reviewed, and the slit function and splitting correc- 
tions are put on a sounder mathematical basis without changing the re- 
sults. Additional effects and assumptions are considered. 


Line Shapes 

In most charge-exchange collisions involving electron exchange, the 
velocity of the projectile is little affected by the collision (ref. lU). 
Consequently , for a charge-exchange process of the type 

A^ f B > A* + s'" (V) 

where the proJectl_e A* is a positive atomic ion, the target B is a 
neutral atom or molecule. A* is an electronically excited charge- 
exchange neutral atom with the same nucleus as A"*" , and B'*' is B with 
one electron removed, a measurement of the A"*" temperature can bf ob- 
tained by measuring the line width (or second moment) of an A* line 
provided the A'*’ temperature is much greater than the temperature of B. 
The second moment of A** is obtained theoretically by studying the 
Doppler-broadening of the charge-exchange component of one of its spec- 
tral lines. The Doppler-broadening appears in the line shape function. 

The line shape function F is proportional to the emission intensity 
at a specific wavelength for an optically thin line. A typical line shape 
for the charge-exchange component of is the wide bell-shaped curve 

including the two faired segments in figure 5. For analytic purposes it 
is convenient to normalize the line shape so that 
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F(A) dA - 1 (13) 

where A Is the displacement from the line center In terms of wavelength. 


Convolution of Line Shapes and Slit Functiwa 

Ihe apparent line shape F observed at the exit slit of a mono- 

cnromator as the wavelength setting of the monochromator is changed is 
actually the convolution of three shapes: (l) the line shape F^ of a 

fine structure component (they will be assumed to all have the same 
shape); (2) the line shape F^ of the fine structure at zero tempera- 
ture (actually a series of sharp peaks); euid (3) the slit function S 
of the monochromator. The line shape F^ is what is needed to deter- 
mine the temperature. 

The line shape F^ is given by 

Fj^(A) - I^6(A - A^^) (lU) 

where is the Integrated Intensity (the spectral intensity integrated 

over all wavelengths) of each component, 5 la the Dirac delta function, 

A^^ is the displacement of the center of the i^^ component from the center 

of the line, and the sum is over the k components. Equation (1^) in- 
cludes both the usual fine structure cmd magnetic splitting. The shape 
F^ satisfies equation (13) if the I^'s are normalized by 


I. - 1 (15) 




The monochromator slit function S accounts for line broadening 
that is due to the finite width of the entrance and exit slits , optical 
aberrations, and unwanted diffraction. The slit function is proportional 
to the response of the monochromator when set at one wavelength but il- 
luminated by a different wavelength. Two typical slit functions are 
shown in ligure 12. The slit function is normalized by 



dx » 1 


(16) 


A 


10 


Tht true lint shape is the convolution of and F^ (ref. 15). 

Fj.(A) - ^ F^{h - A')F^(A') dA' (17) 

The apparent line shape is the convolution of F^ and S (ref. l6). 


F^(A) • ^ S(A - A")F^(A") dA" (l8) 

Substituting equation (17) into equation (l8) and reversing the order of 
integration gives 


^ F^(A’) ^ S(A - h")F^{h" - A') dA" dA' 


(19) 


Making the substitutions A' 5 A - A"' and A - A" = A"" 

(A - A'") / S(A"")F2(A"' - A"") dA"" dA" ' 




( 20 ) 


But the inner integral is simply the convolution of S and F^, while 

the outer integral is the convolution of the inner integral and F^^. 

Hence F^(A) is independent of the order of convolution. This Justifies 

the order in which the second moment corrections are made in reference 2 
and in appendix C of this report. 


Second Moment in the Limit of Zero Ion-Cyclotron Radius 

In the determination of ion temperature from F^, it is desirable 

to correct for monochromator slit function, fine structure and magnetic 
splitting, and the superposition of charge-exchange spectra from several 
reactions. These corrections are greatly facilitated if the Doppler 
broadening is measured in terms of the second moment insteaui of the half 
width of F. The second moment is defined by 




A^F dA 


( 21 ) 
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which is the average value of the square of the wavelength deviation of 
the line shape from its center. 

Ho loss cones or drift . - It is shown in reference 2 that in the 
limit of zero ion-cyclotron radius and assuming spatial homogeneity and 
no loss cone or azimuthal drift that the second moment is 



( 22 ) 


where the symbols are defined in appendix A, and is the second 

moment of F-. Cross sections for H and H. are shown in figure 13. 
d a D 

Results from equation (22) are tabulated in reference 2 for Hg, 

He 587.6, and He 667.6 lines and are plotted in figure 6 of this report 

for H , H., D , D. , He 507.6, and He 667.6 lines. In the case of the 
a B a 6 

Balmer lines, the target vas assumed to be H2 or D^. 

Loss cones but r.o drift . - If there is no azimuthal drift but ions 
are lost through loss cones in velocity space, the second moment can be 
derived from the relation 



(23) 


by assuming a Maxwell-Bolt zmann velocity distribution outside the loss 
cones. The result is 



(2U) 


where R is the mirror ratio.. This is simply (1 ♦ 1/2R) times equa- 
tion (22). Because R > 1, this means that (a^i. will be increased. 

Hence, if is determined f rom will be smaller if equa- 

tion (2U) is used instead of equation (22). 
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Drift but no loss coneg . - It ia aasuned that the plaama ii spatially 
hofflogeneous , that there is a Maxvell>Boltzmann distribution of ion veloc- 
ities about a uniform azimuthal drift, and that the neutrals remain sta- 
tionary. If the plasma is viewed across a diameter, and the ion cyclotron 
radius is zero, the azimuthal drift can be replaced by a drift in the 

direction (fig. lU), where the line of sight and i>agnetic axis are in the 

V and V directions, respectively. The case of drift but no loss cones 
X z 

can be derived by the same general method as in appendix B with the result 




o,(E)e (cosh Q - 


s 


-E/kT. . , 

o,(E)e ^E dE 

4 


dE 


(25) 


Herce, in general, drift affects the second moment through its influence 
on E, o„, and Q. The exception to this is the special case where 

o„ % v~^ so that the charge-exchange neutrals have a Maxwell-Bolt zmann 

velocity distribution like the ions. In this case if observations are 
made perpendicular to the drift velocity (as was assumed in deriving 
(25) )• the second moment is independent of drift in the zero-ion- 
cyclotron radius approximation used here. However, in practical cases 

o„ does not vary as v”^. In fact, for o„ Increases with increas- 
ing V (see fig. 13). The second mcment for is shown in figure 15 

for no loss cones and various drift velocities. For this cross section 
(fig. 13) drift tends to reduce the second moment. If the drift becomes 
much larger than the ion thermal velocities, the v in o„(v)v is 

almost independent of the component of v along the line of sight, so 
o„(v)v may be regarded as a constant. Hence the second moment approaches 

the value for a Maxwell-Bolt zmann distribution at the ion temperature. 

Loss cones and drift . - If there is azimuthal drift, the loss cones 
are reduced in size and are no longer conical. The same basic assump- 
tions as in the preceding section are made. This case is treated in 
appendix B. Based on floating probe and magnetic flux density measure- 
ments in SUMMA, a typical drift velocity is UOO 000 m/sec. t’or this 
drift velocity, the effect of mirror ratio on shown in fig- 

ure l6. Here the value 1.66 of the mirror ratio is the average of the 
mirror ratios of HIP-1 and SUMMA. Just as for the zero drift case, 
the second moment is increased by the presence of loss cones. The 
second moment for H^ is shown in figure 17 for a mirror ratio of 



13 


1.68 and varioua drift velocities. The drift tends to reduce the second 
oooent even more than in figure 15. If the drift is large enough, it 
mcdies the loss cones negligible, and the second moment approaches the 
t4axwell>boltzmann limit Just as in the previous case. 


Correcvions to Second Moment for Slit Function 

A correction c to the second moment due to monochromator slit func- 
tion is derived in appendix D of reference 2 for a Gaussian slit function 
(see fig. 12) and is 


e - - - -0.1803 (26) 

where W is the full half width of the slit function. For many mono- 
chromators this Ic a good approximation for narrow entrance and exit slits 
of the same width. 

To obtain a good signal-to-nolse ;j-1o, it is frequently necessary to 
use wide slits. For the ideal case of wide slits of equal width, a tri- 
angular slit function is produced (see fig. 12). In appendix C it is shown 
that the second moment correction for a triangular :<lit function is 

e ■ - ^ ■ -O.I 667 ( 27 ) 

Assuming the slit function correction is considerably smaller than the 
second moment, an average of equations ( 26 ) and (27) namely -0.1735 
should be accurate enough for Gaussian or triarguleur slit functions or 
slit functions intermediate between the two. 


Second Moment for the General Caise 

In reaction (V) the target B in the case of hydrogen could be 

either H or Hg. The projectile A* could be either or Hg. 

In addition to reaction (v), twe-step processes (ref. 2) may be important 
if the plasma is large enough or dense enough. In general, the second 
moments for different processes will be different. It is shown in 
appendix D that for n charge-exchange processes with superimposed 
charge-exchange spectra the second moment total line shape 


L 
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('•')3.t ■ B'j 3 • *oj) 

J-1 

where Ij is the integrated intensity for process J. 

corrections are made for splitting (app<“ndix C of ref. 
mator slit function 

E ‘j * ''oj) • ■ E 'i"oi - 

J-l i-l 

where is the second noment observed at the exit slit of the mono- 
chromator as wavelength is scanned. Equation (29) can be solved for 

if a table of the left hand side versus is calculated and inter- 

polated. 


( 28 ) 

Consequently, if 
2 ) and monochro- 


Effect of Finite Ion-Cyclotron Radius on Half Width 

In the above derivations it was assumed that the ion-cyclotron 
radius was zero, although equation (29) is applicable even for finite ion 
cyclotron radius. Recently the problem of finite ion cyclotron radius 
has been treated (ref. 3) using the guiding- center model. It was assumed 
that (l) there was *d loss cone, (2) the ions had a diffusion density 
profile in the radi.al direction, (3) the neutral concentration was con- 
stant, (1*) the drift velocity of ions was e * 2/B^, (5) the Ion tempera- 
ture was a constant, and (6) the charge-exchange spectra were due to 

+ H2 H* + H* (VI ) 

The combination of finite lon-cyclotxon radius with azimuthal drift 
was found to increase the half width and hence the second moment (the 
second moment is roughly proportional to the square of the half width). 
This is Just the opposite of the effect of drift in the case of zero ion 
cyclotron radius (figs. 15 and 17). The finite radius effect can be 
understood from figure I8, which assumes no delay in emission by the 
cheurge-excheuige neutral. The ion in the line of sight has the usual 
cyclotron velocity Vg. However, its guiding center is located above it 
and to the right and has a drift velocity v^. When v^ is added to 

V , the resultant velocity v of the ion has a larger component along 
c r 

the line of sight than does v^ alone. Hence the spectral line is 
broadened. The broadening effect of drift due to finite cyclotron radius 
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waa found dominant over the ov effect for the relatively moderate ranges 
of drift to random energy reported in reference 0. 


Effect of Cascade on Half Width 

In reference 8 the effect of delayed cascading from higher levels and 

of delayed emission of were also Investigated. These had little effect 

on the half width of the emission observed across a diameter, but caused 
large Increases in half width for lines of sight that did not pass through 
the center of the plasma. Unfortunately, insufficient cross sections are 
available for hydrogen to evaluate cascade and delayed emission properly. 


Inversion 

It would be desirable to observe emission for a number of values of 
y (fig. U) and then invert the spectra to obtain as a function of 

radius. Unfortunately, an Abel inversion is not applicable to the spectral 
intensity because the spectral intensity is anisotropic. However, the 
Integrated intensity (area under the spectral line profile) is isotropic 
so an Abel inversion is applicable to it. Nevertheless, there are ob- 
stacles to obtaining as a function of radius: (l) cascading and 

delayed emission in hydrogen must be better evaluated, (2) it is not 
known with certainty wnat formula is a satisfactory approximation to the 
Ion azimuthal drift velocity, and (3) numerical methods must be developed 
for the inversion itself. It appears that the second moment is better 
suited to this problem than the half width. 


CONCLUDING REMARKS 

The first purpose of the study was to obtain scaling relations for 
ion temperature by a statistical method. The ion temperatures were meas- 
ured in nonequilibrium hot-ion plasmas by observing second moments of 

the charge-exchange neutral components of atomic linej. These moments 
were then corrected for three effects and compared with moments calculated 
for vai’iouB ion temperatures. In SUMMA, midplane ion temperatures ranged 
from 1£9 to 05^ eV in hydrogen and from 666 to 1*927 eV in helium. In 
HIP-1, midplane ion teraperatvires were from 170 to 607 eV in hydrogen, 
from 260 to 650 eV in deuterium, and from 597 to 900 eV In helium. The 
scaling relations were 

-V j0.52±0.10yl.61*±0.05g-l. 3010.22 g+j 
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-V j-0*li^0*2\l*36±0.05g-0.29l0.22 
^ j0.33±O.OUyl.50t0.10j^-0. 1*910.10 

-v i0*32l0.05g-0.30l0.12 (gjp.j^ 

No scaling relation was obtained for helium in HIP-1 due to insuffi- 
cient data. After the power law relation has been established, the 
statistical method used has the advauitage that it is not necessary to hold 
two variables constant while varying a third when taking the data. 

The second purpose of the study was to advance the theory of Doppler 
broadening of charge-exchange neutral spectra in rotating plasmas in 
mirror devices. It was shown that the order of making corrections is im- 
material. Loss cones were found to increase the second moment. In the 
limit of zero ion-cyclotron radius, azimuthal drift decreetsed the second 
moment when the plasma was viewed across a diameter. 'Riis effect is 
opposite to euid was obscured by the influence of finite cyclotron radius 
in the range of the ratio of drift to random energy reported previously. 

To obtain more accurate ion temperatures, drift will have to be included 
in the future. 

Hie correction for a triangular slit function was found to be about 
the same as for a Gaussleui slit function of the same width. 

The use of the second moment Instead of half width should facilitate 
obtaining the ion temperature where several kinds of charge-exchange 
processes produce a single charge-exchange line component. 


APPENDIX A 


UST OF SYMBOLS 
diumy variable 

midplane magnetic flux density 
coefficients of net regx'esslon 
velocity of light 

Integral In denominator of for loss cone and drift 

first term In D 

second term In D 

total kinetic energy of Ions 

electric field 

line shape function 

line shape function observed at exit slit 

F for A - X < W 

e mo 

F for X - X > W 

e mo 

true line shape function 

line shape function for fine structure at zero temper- 
ature 

line shape fiinctlon for fine structure component 
Wc(M/2kT)^^^/X 

o 

v.(2M)^^^/kT 

a 

integrated intensity (spectral intensity integrated over 
all wavelengths ) 

Boltzmann's constant 

mass of Ion 

integral In numerator of for loss cone and drift 

first term in N 
second term in N 

number of charge-exchange processes 
v^(M/2)^^^ 


mirror ratio 
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r 



U 

V 


V 

V 

c 



V 

p 


V ,v ,v 
X* y* z 


w 

X 


Y 

y 


z 

Y 

6 


e 

e 


cl 

*c2 


• ’ • ,r'J 


m 



A„c(M/2kT)^^^/A„ 
n o 

slit function of monochromator 
standard error of estimate of 

standard error of estimate of z 

midplane Ion temperature 

(R - 1)P/E^^^R 

electrode voltage 

velocity 

cyclotron velocity 
drift velocity 

component of velocity peirallel to 5 
resultauit velocity 
components of v 

full half width at half intensity of S in wavelength units 
E/kT 

difference between wavelength of light and monochromator 
setting, X - X 

n 

(R - 2(R - 1)P^/E)^^^/R 

disccmce from line of sight to centerline of test section 
Uee fig. U) 

Mc^A^/2X^ 

o 

Dirac delta function 

correction to second moment due to monochromator slit function 
polar angle in velocity space 
9 for top of loss "cone'* surface 
0 for bottom of loss "cone" surface 

displacement from center of charge-exchange component of 
spectral line, X - X^ 

difference between wavelength set on monochromator and center 

of line, X - X 
m o 

displacement of center of component frcm center of line 

apparent second moment of charge-exchange component of spectral 
line 
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corrected for monochromator elit function, magnetic splitting , 
and fine structure 
X wavelength of light 

X_ monochromator setting 

IB 

X^ wavelength of center of charge-exchange component of spectral line 

0 ^ charge-e/.change cross section for excitation to electronically 

excic'd state 

^ azimuthal angle in velocity space 

♦ for loss "cone" surface 
Subscripts : 
f fast 

i fine structure component index 

J process index 

k number of fine structure components 

n nximber of processes 

t total 


Special conventions : 



average 


I 
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APPENDIX B 

DHIIVATION OF SECOND MOMEWT FOR DOPPLER-BROADENED CHARGE-EXCHANGE 
NEUTRAL SPECTRA AT THE MAGNETIC MIDPLANE 
WITH DRIFT AND LOSS CONES 

When azimuthal drift ii present in a magnetic mirror apparatus, the 
size of the loss cones is reduced (ref. 17) and they are no longer conical. 
Since the population of the loss cones is (generally assumed to be negligible, 
and since the neutrals are assumed not to drift azimuthally, the second 
moment of charge-exchange spectra is not the sane as for no drift and no 
loss cones even if the line of sight is across a diameter, as ve will assume. 
We will also assume that the ion-cyclotron radius is zero (or, at least, is 
negligible). 

From reference 17 the loss "cone" surfaces are defined at the magnetic 
midplane by 


- v^(R - 1) + (Bl) 

pc d n 

where is the component of ion and charge-exchange neutral velocity 

parallel to the magnetic field, v^ is the ion cyclotron velocity, and 

v^ is the ion azimuthad drift velocity (the charge-exchange neutral has 

the same velocities as the ion at the time of reaction but v becomes 

c 

constauit in direction). We adopt polar coordinates as shown in figure lU. 
Letting 


^ ^ ^dV2 


(B2) 


and 


E « j Mv^ 

equation (Bl) can be solved for the value of ^ on the loss "cone." 

\J2 


(B3) 


s sin 
c 


-1 


1 - 




(R - l)(co8^e Yj 2 ® ” R 

£ 


f 1 pf_\ 
R E / 


1 - cos e 


(BL) 


21 


There are obvloualy limiting values of 6 outside of which there are 
no loss cones. These are found by setting ♦g ■ 0 in equation (BU) so 


where 


and 


cos 0 , ■ U 
cl 

Y 

(B5) 

cos 0^2 ■ U * 

Y 

(B6) 


U 


iR-l2P 


(B7) 


Y 


R - 2(R - 1) 



R 


(B8) 


There is a limiting value of E (or v) below which there are no 
loss "cones." This is found by setting 0^^^ ■ 0^^ in equations (B5) and 

(B6) with the result 


2(R - 1)P^ 
R 

The second moment is given by 



(B9) 


(BIO) 


where the average is over the charge-exchange neutral velocity distribu- 
tion function, the unnormalized version of which is assumed to be (in 
Cartesian coordinates and for plasma along the line of sight and to the 
right of the vertical centerplane (fig. i+)) 


.(v) 


ve 


-[M/2kT][v^(Vy-v^ 




(Bll) 


which is a Majcwell-Boltzmann distribution with a drift correction and a 
a„v correction. From equations (BIO) and (Bll) 
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(B12) 


where 


H 3 


v^(2M) 


1/2 




(B13) 


If we call the triple integrals in the numerator and denominator of equa- 
tion (B12) N and D, respectively, then 


2A 

/,2\ o N 


(Bli*) 


By symmetry it is only necessary to integrate N and D over the range 

0 < 6 < TT 


0 1 ♦ < ^ 


0 < E < • 

excluding the loss cones. This gives 

itN, N_ 

« • -i* TT 

2H 


(B15) 


where 


"‘■‘I 


Eo,,e"^^^'^(cosh Q 


^osh Q - dE 


(B16) 
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t\ [i |.‘‘''r2(u" * y" - 1) ,l„h Z 

♦ M (u - i) (cosh Z - Siiip)] • i (co.h (» - 

(l - 008^6 - 2^ 

c 

♦ iln (2<>g)] d (cos e)j dE 

where 

Q 5 
Z 5 

Similarly 


+ e 


-E/kT 


'-'eo.e^j 



where 




-E/kT 

Q * 


sinh Q d£ 



(BIT) 


(B18) 

(B19) 


(B20) 


(B21) 


(B22) 



2k 


Hence (a) is calculated fron equations (B2), {Bh) to (B9)i and (B13) 

to (B22). For small arguments, expressions of the form 

cosh A • slrh A/A should be expanded In power series to avoid loss of 

accuracy. 

The above derivation Is for the plasma to the right of the vertical 
center plane. However, along a diameter the plasma to the left of the 
vertical center plane produces the same symmetrical unshifted line shape 
and second moment as the pi earns to the right. 
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APPilHDlX C 

DERIVATION OF SECOND MOMENT CORRECTION DUE TO 

TRIANGULAR SLIT FUNCTION OF MONOCHROMATOR 

Tht apparent second moment should be corrected for the slit function 
of the monochromator. We assume a triangular slit function S (fig. 12) 
for a monochromator set at X . The slit function is u/mmetxic about X 

Q D 

and is given by 

A - X ♦ W 

S . 2 (X - W < X < X ) (Cl) 

,2 u _ — a 


X - X ♦ w 

m 


(X ♦ W > X > X ) 
m — — m 


We assume the line shape F before pausing through the monochromator is 
the shape for a constant charge-exchange cross section (since the order 
of convolutions does not matter, as proved in the section IMPROVED THEORY 
OF LINE SHAPES AND SECOND MOJENTS, fine structure and/or magnetic splitting 
can be hypothetically removed before passing throu^ the monochromator). 
From equation (B7) of reference 2 this is 


.,1/2 

F(X - X ) - 

X (2kT) 
o 




2 -w . 
we '^w 


(yAt)^ 


where ■ E/kT. From equation (D6) of reference 2 the apparent line 
shape is 


e m o 


F(X - X^)S(X^ - X) dX 
o m 


Substituting equations (Cl) to (C3) into (CU) gives 


F.(X - X ) - 

° VTX (2kT) 
o 


O m O 

172 / 

o. -w vJ(yAt)- 


w e dw(X - X W) dX 

D1 


>X +W rN' 


p 

we dw(X - X ♦ W) dX 
zn 


(y/ht)- 
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Thil will be evaluated only for the caae because the case 

X < X will be synsnetrle. For X > X there are two subcases as 

m — o m ~ o 

sho'.n In figure 19. For cases (a) and (b) we designate the apparent 
line ch .J by F and F g* respectively. 


„l/2 

F , (a - X ) - 2- 

° W^X (2kT) 
o 


_vJx -w vJ-{ 




0 

we dw(X - X ♦ W) dX 

IB 


rs. 




w^e dw(X - X ♦ W) dX 
n 


. r‘"*“ r 

Jx vJ(Y/kT)^^^ 

n 

“ »o(2“T) U X -w J(r 


'.2‘\ - ‘o> 


p ! 

we dw(X - X ♦ W) dX 

ID 


(yAt)^^^ 


w^e"'' dw(X - X ♦ kWX 
m 


>(y/kT)^^^ 


? — w 

we dw(X - X ♦ 

IB 
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and 



So equation (C8) becomes 



This can be integrated with difficulty to give 



But this is simply 

{<} ■ * 4 

wnere (a^)^ is the second moment of F. We assume equation (C13) 
good approximation even if is not exactly constant. 


(CIO) 


(Cll) 


(C12) 


(C13) 
is a 
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APPENDIX D 

DIBIVATION OF SECOND MOMENT FOR DOPPLER-BROADENED CHARGE- EXCHANGE 
NEUTOAL SPECTRA WITH SEVERAL SUPERIMPOSED PROCESSES 

We assume that the line shapes Fj of several charge-exchange 

processes J ■ l,2,...,n are superimposed. We do not assume that the 
F. have the same shape, same uecond moments, or same centers. It is 

J , 2 \ 

desired to find the second manent (A V of the line shape resulting 
from superposition. 

The Intensity of component J resulting from process J at 
wavelength A is proportional . integrated 

Intensity of component J normalized by 




The wavelength A of the center of the line shape for the superposition 
is ° 


' Up j»i J J 

J r\<» n 

. S"*' 


We define the center of profile J by 


^0 


so equation (D2) becomes 


A . I.A , 

o Zy J oJ 
J-i 
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Obviously 



If we define the second moment of component J by 



(X - dX 

oj J 


and define 



then equation (D^) becomes 




J-1 


(D5) 


(d6) 


(D7) 


(d8) 
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TABLE I, - RANGE OF TEST CONDITIONS AND ION TEMPERATURES 


Apparatus 

Gas 

Midplane -aagnetic 
flux density, 

B, 

tesla 

Electrode 
voltage , 
V. 

KV 

Total cur- 
rent, 

I, 

A 

Pressure of 
neutrals , 
torr 

Midplane ion 
temperature, 

^i‘ 

eV 

SUMMA 

Ho 

2.07 to 2.U9 

7.1 to 19.6 

0.60 to 1.10 

8.8»10"^ to l.l-lO""* 

169 to 85^ 


2 

He 

2.07 to 2.U6 

U.9 to 19.9 

0.82 to 1.06 

-J* -1* 

1.3*10 to 2.0-10 

666 to U927 

HIP-1 

“2 

0.75 to 1.12 

12.0 to 22. U 

0.U9 to 1.92 

5.U«10"^ tc 1.3-10“ 

170 to 607 


Do 

0.60 to l.lU 

22.0 

0.21 to 1.08 

6.9-10 ’ to 1.0-10 

268 to 650 

1 

2 

He 

0.81 to 1.10 

17.0 to 22.0 

0.38 to 1.28 

-U -U 

l.U-10 to 1.7-10 

597 to 980 


w 




TABLE II. - STATISTICAL SUJ«ABY OF SCALIMI RELATIOBS 


Apparatus 

Gas 

Sample 

size 

Independent 

variable 

Sample partial 
correlation 
coefficient 

Level of 
signifi- 
cance 

Correlation 

significant? 

Standard devia- 
tion of coef- 
ficient of net 
regression 

Standard 
error of 
estimate 
of log 
T 

S 

z 

SUI-IMA 

H 

2h 

Overall 





0.0298 




I 

0.760 

<0.001 

Highly significant 

0.10 





V 

0.989 

<0.001 

Highly significant 

0.05 





3 

- 0.801 

<0.001 

Highly significant 

0.22 


JUMMA 

He 

38 

Overall 





O.OU 53 




I 

-'^.076 

0.663 

Hot significant 

0.2i* 





V 

0.976 

<0.001 

Highly significant 

0.05 





B 

- 0.216 

0.205 

Not significant 

0.22 


HIP-1 


88 

Overall 





0.0J*85 




1 

0.687 

<0.001 

Highly significant 

O.OU 





V 

0.855 

•0.001 

Highly significant 

0.10 





B 

- 0 . U 72 

<0.001 

Highly significant 

0.10 


HIP-1 

Dp 

U2 

Overall 





0.0593 




1 

Hot a^-ailable 

<0.001 

Highly significant 

0.05 





B 


0.016 

Significant 

0.12 
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Figure 1. - Test section and magnetic field configuration for SUM^TA. 
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Figure S. - Typical line exhibiting wide and narrow components. 
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